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SYIOPSIS 


In this work th.e problem of deterinning tie maximal 
arc-dis joint axio node -disjoint flow in a multicommodity 
network is considered. This is the restriction of the 
well known multicommodity maximal flow problem with the 
added constraints that each arc allows flow of at most 
one commodity ( arc-dis joint flow), or each node allows 
flow of at most one commodity (node-dis joint flow). 

Both L'he pi*oblems for two-commodi cy airected 
necworbsare formulated as integer linear programraing 
problems and methods are suggested to solve them. Branch 
and bound methods are developed to obtain optimal solu- 
tions, and heuiistics to obtain approximate solutions. 
Computer programs are designed for each of the method 
suggested and tested over a large number of randomly 
generated networks. Computational performance of the 
programs on iii;C“1090 time-sharing computer system is 
presented. 

Arc-dis joint and node-dis join t floors in multi- 
commodity networks are also considered. Two heuristics 
are developed for each of the problem. Computational 
performance of the heuristics over randomly generated 
problems is also presented. 



CHAPTER I 


INTRODUCTION 


THe aim of an operations researcher lias always been 
to tackle bhe unbackled problems, to solve the unsolved 
problems ana to propose better solution procedures for the 
solved ones. This work is also consistenx with this aim. 

I have tackled unGaclcled problems, proposed methods to solve 
them and givoi hints for their improvemenbs and extensions. 

On account of its wide applicability, network flows 
are considered an important branch of scuay in the field of 
Operations Research. Problems in the study of airline 
systems, railroad systems, shipping syscems and in general 
all transportation planning problems are network problems 
on appropriately defined networks. Problems an natural 
gas, crude oil or other fluid flows are problems in the 
sppropriate pipeline networks. Economic models can also be 
treated as network models by representing factorieo, ware- 
houses and markets as points and by treating highways, 
railroads and wameiways, and other transportation channels 
as lines. Thus nexwork algorithms find a vast number of 
applications in co li unication, transportation and distribu- 
tion and in solving various practical problems that can be 
modelled as network flow problems. 



If a network has many source nodes and many sink 
nodes and distincx commodities flowing be'cireen each source- 
sink pair, then it is a multicommodity netirork and flow 
problems over it are referred to multicommodity network 
flOT,r problems, due to technological reasons if arcs allow 
only one commooitAr to flow over it, then any feasible flow 
T/jill be arc-dis joiiit , i.e., commodities will flow over 
disjoint (mutuallj'’ exclusive) sets of arcs. Similarly, 
node-dis joint flow can be defined as flow- in winch commo- 
dities flow over disjoint set of nodes. 

The problem o£ determining maximal a,rc-dis joint and 
node-dis joint flow in a multicommodity network is considered 
in this work. The problem may be looked upon as to sub- 
divide the set of arcs (or nodes) into as maxiy subsets as 
the number of commodities. With each comiodity flowing over 
its subset of arcs (or nodes), ib is required to find that 
subdivision which enas up with maximum sum of flows. 
Evidently it is a combinatorial optimizauion problem Tcith 
diflicult solvable properties. One has bo device special 
solution techniques to solve them. 

Chapter II presents a state-of-arx survey of results 
and algorithms for multicommodity flow problems. This 
survey is supplemenxed by an extensive list of refei*ences. 

Chapter III considers edge-dis joint maximal flow 
problem in a two-commodity directed ne-Lwork, Ib describes 
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integer linear programming formulation ox the problem, 
retails of a braiich and bound algorithm developea to solve 
bins problem optimally and of two heuristics developed to 
solve large sized problems are described. Computational 
performance of all the algorithms descx-ibed in this chapter 
is also presexXi-ed, 

Chapter IV considers maicimal node~dis joint flow in a 
two-commodity direcced network. It descx’ibes the mathema- 
tical programming^ formulation of the problem; details of a 
branch and bound algorithm to solve the problem optimally 
and of two heurisiics to solve lai'ge sized problems. In 
essence edge -dis joint flow problems and node-dis joint flow 
problems are sixiiilar problems, hence solution procedures for 
them are along similar lines. 

Extensions of the heuristics developed for two-commo- 
dity neiworks to multicommodity networks is considered in 
Chapter V. Compucor programs for them ore also ifjritten 
end computatioxial performances of the programs is discussed. 
Lasuly, in Chapter VI avenues for further research are 
mentioned. 



CHAP'iEE II 


A REVIE¥ OE IFJLTI COmOEITY HETWORE I'10¥ PROBLEMS 


'This chapter aims at a compreliensive survey oL the 
literature cLealiii£ with the multicoimodi fcy flow problems. 

Ihe survey may loj,ically he spread into three sections. 

Section 1 discusses the properties of mufL'cicommoctity net- 
work flow problems. Section 2 presents a brief survey of 
special results. General multicomuiodi cy network flow problems 
and their solution procedures are discussed in Section 5. 

2.1 G-ERERAL RSLARKS 

Multicommodicy flow problems arise naturally in 
modelling wherever commodities} vehicles or messages are to 
be shipped or tranbnaicted from certain nodes of an under- 
lying network xo some others. Considers cions to simple 
examples make it clear that the multicommodity flow problems 
are considerably more complex xhan the single commodity one. 
I'he nice combmacorial features of che single commodity case 
are lost in its generalization. The coiic tramt matrix is 
no longer unimodular hence optimal solutions are generally 
non-integer, flie max-flow, min-cut theorem true for single 
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commodity ne CT7orlcs is false for mulxicoinnoaity netvorks and 
no simple minded modification of the leloelling process seem 
CO work. Ho we wo 1 j che constraint matrix lias a special 
structure of heiiip block angular which ic exploited to 
produce efficient: algorithms. 

2.2 SPECIil EHSblfS 

Most of the special results for mulfcicommodixy network 
flow problems arc restricted to special networks or netwox'lcs 
with at most two sources or two sinks or cwo or tliree 
commodities, fhese results provide litble aid in actacking 
the general problems, hovevex, their study is advantageous 
before proceeding to the general problems. 

2.2.1 fw o Co mno di t y Plow Pro blem s 

l‘¥ 0 -commo( itj' flow in undirectea no cworks enjoys many 
special properties. Max-flow, mi/i-cut thcoci^m is valid for it g 
hence optimal solutions are integer. Algorithms for 
constructing ma^^imal two- commodity/ flo'vr axe efiicient and 
may be found in [16,18,30]. 

Two-comraob ity flow in directed networks does not 
possess the properties of two-commodity flow in undiiected 
networks. OptiiiiaA solutions are generaJ.ly non-integer. It 
IS harder to solve tins class of problei'is, in lacCj, as 
difficult as linear programs. 
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2*2.2 Comaon S ia k and Common. Source Mu lti comodity Plow 
Problems ' * ’ 

A mixlticommociity network, which has exactly one sink, 
i.e., all commod-i Lies xerininate at one point, can be solved, 
very efficient!-^ due to iis special structure. Max -flow, 
min-cut theorem is valid and hence opbimal values are 
integer. Different commodities may bo created as a single 
commoaity while solving the problem, and later flow values 
of each commodi cy may be identified by an arc-chain 
decomposition. Algorithm for finding maximal flow is 
presented in [29j. Minimum cost flow problems for these 
networks can also be solved efficiently employing the same 
ideas. 

Common source multicommodity flow problems can be 
transformed to common sink multicommodicy flow- problems by 
a Simple transformation and then all efficient algorithms 
available for com on sink problems may be applied to it. 

2.2.3 hli s ce 1 1 one Q u s Re su it s 

Rothfarb and Prish [28] present a max- flow, min-cut 
cheorem and an algorithm for a special six-node, three 
commodity flow problem on a complete graph. Eleitman [25] 
presents an algorithm to obtain maximal flow when each 
node of the graph is a sink for all bub one commodity. 
Recently Evans [7] has described suflicienb condition by 
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which a certain class of multi commodity minimum cost 
network flow prollem can be transformed and solved as an 
equivalent single commodity flow problem. 

2.3 G-EKERAL MULJI COMMODITY FLOW PROBLEKS 

This sec'Jion discusses some important multicommodity 
flow problems and xheir solubion procedures. 

2.3.1 Multi commoa it y Network Realizat i on Problem 

The problems concerning that prescribed flow values 
of commodities can be realized in a given network or not 
come under tbas heading. These probloms are encountered 
quite occasionally in practise and hence researchers have 
attempted to taxhle it. One solution technique for tins 
problem is presented in [16]. It is ba&ically a revised 
Simplex technique and it may solve moaerately sized problems. 

2.3.2 Networ k S ynthesis Problem 

In this problem we seek an optimal network configura- 
tion sacisfying a given set of requirements. These problems 
may be viewed as specializations to ne ttrork desigxi problems. 
In general, ire expect them to be more difficult sanice each 
tentative configuration may have to be analysed lo find the 
optimal configuration. A survey of li cera'cure available 
for this class of problems is given in [36], 
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2*3.3 Multi CO fi ffliodity Maximal Plow problem 

jjetermiiiijag maximum sum of flo\7 of all commodibies in 
a mulGicommodi!:3;- network vrliicli satisfies tlie capacity 
constraints a.n(l node conservation consbraints belongs to 
tins class of problems. Pord and Fulkerson [9] xj'ere 
probably the first to formulate this problem as a linear 
progrs'-i and propose solution techniques. They proposed the 
arC“Chain formulation of the problem. I'he problem is solved 
by revisectsiiiiple.L method which requires laaowledge of the 
basis only. I'hc variable to enter the basis is determined 
by solving sho'>’test pa.th problem for each coramodity. 
Unfortunately no compubat lonal results are available for 
this technique. 

I'he same problem has been formulaGed as a node-arc 
formulabion in [19] and solvea by the xiiethoa simlar to 
that of arc-chain formulation. Grinola [14] has suggesbed 
a resoui’ce direcbive solution technique which is reported' 
to have slow convergence proper bies. 

2.3.4 Mul-tic on loa ity Minimum Cos b Fl oir P roblem 

I'he multicommodity minimal cosb flow problei^a is to 
determine the Eniiimum cost mulbi commodity floir through the 
network that meebs the demand of each coinmoditj/ j subject 
to (l) supply res Grict ions ; (2) arc-capacity restrictions ^ 
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and (3) flow conservation at each translnpment node. This 
problem is a linear program and may be solved by the simplex 
method or one of its variant. Howevex'', real world problems 
are frequently of such size xhat direct application of 
simplex 13 prolnbitive. There are bhi''ee basic approaches 
that have been proposed for developing specialized tech- 
niques for these problems. They are price-directive 
decompositions resource -directive decomposition and parti- 
tioning methods. 

A price-direct ive decomposition pi’ocedure directs the 
coordination between a master program and each of sevei-al 
subprograms. The objective is to obtain a se'c of prices 
such that uhe combined solution for all subproblems yields 
an optiiiial for the original problem. The master program is 
solved by the revised simplex method xrich the subproblems 
used to test for opcimality and select candidaxes for 
entering the basis of the master. The subproblems are 
single commodi cy minimum cost flow problems and can be 
solved efficiently. These procedures fox- the original 
problem and igs generalizations may found in [ 3? 31 , 34, 35 ] • 

A reso’Tce dii'’ective decomposition procedure; -when 
applied to a multicommodity network floxr problem haviiXg 
q commodities ivO to distribute the ai'^c capa.city among the 
indivisual commoaities in such a way that solving q sub- 
problems yields an optimal flow for the coupled problem. 
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A't- eacli iteration an allocation is made and q single com- 
modity flow problems are solved, flie sum of the capacities 
allocated to an axc over all commodities is equal to the 
arc capacity all oca cion is developed. Robacker [27] was 
the first one to propose these ideas, however, he did not 
inuicate how xlm c approach could be implemented. G-eoffrion 
[10], Ilennington and Shalaby [24] and Eeld, Rolf and 
Crowder [15j have presented techniques for implementing 
this approach. 

PartitiOiiin^ xechniques are specialization of the 
simplex method irhere constraints are partirioned into two 
sets known as current constraints and secondary constraints. 
I'he secondary constraints are relaxed and objective function 
IS minimized Trliile only current constraixi'cs are active. 

If the optimal solution thus obtains a sa.tisfies the secondary 
constraints, ic is an optimal one for che original problem. 
Otherwise those secondary constraints, which are /lot satis- 
fied by the cin’ronc oprimal solution, are added to the set 
of current consti-^ain ts . The current constraints ^rhich are 
not binding at che optimal solution are relaxed and steps 
are repeated. 

The motivation to apply this technique is the fact 
that only a few capacity constraints are binding at opti- 
mality. The paper of Graves and Mcbride [12] deal with 



primal technique while the work of Gn^oriadis and White 
[ 13 ] IS a dual technique. 

It seems ohax field of arc-dis joint and node~dis joint 
flow problems romains uniouched as none of the articles on 
this subject could be founa in literature. 'this is the 
reason that literature survey on dis joint-flow problems 
was not possible. 



CHAPTER III 


MAUI-iAl PIO'/ li. AHC-JJISJOIHT lyO-GOniOuITY NETVORES 


In this chapter we shall discuss arc-disjoint two- 
coramodity m8>’'iiial flow prohlein in a directed network. The 
problem is a v-^riation of the two-commodity maximal flow 
problem in the sense that it is desired that flo^^ of both 
the commodities should be arc-dis joint . These type of 
problems are encouxit ered when the same network (say a pipeline 
system) is used for simultaneous transportation of two com- 
modities and cliG commodities are such that they can not 
flow on the sans arc at same time. 

The analjrsis of these problems is heavily aependent 
on f 0 1 1 owing al g o r 1 1 hms ; 

a) Single coimaotity maximal flow algorithm. 

b) Maximal augmenting paxh algorithm 

c) Branch and bound approach. 

These algorithms are bra efly described in Appendix A. 

This chapter is divided in'co four sections. In sec- 
tion 1 we shall o.iscuss the mathemsticcl programming 
formulation of olie problem as a net^rork flow problem and 
as a 0-1 integer programming problem. A branch and bound 
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method suited for the problem is developed bo get the 
optimal solution of the problem in section 2. As the size 
of the problem uhicli can be solved bj^ branch and bound 
method is limited, tvo heurisiics arc ceveloped anu des- 
cribed in section 3. Compute cional periormancesof the 
branch and bound method and of two heurisxics for I’andomly 
generated nebworlcs is presented in sec cion 4. 

3.1 MAl'HEMAlICii PEOGRAMNG P0 Ri>IULA1I0:j 

The netwonh G(h',A) is a finite set of nodes R and a 
finite set of ordered pair of nodes, called arcs, A. The 
network is sucb that there is at mosx one directed arc 
between any pair of nodes. Bach arc (i,o) has a number 
b attached to it, called capacity, which denotes the 

J- J 

maximum amoimt that can flow on the arc. Two distinct 

11 2 2 

commodities flotr in the network, with s", t and s , t 
being the source and sink node for coimiooity 1 and 2 res- 
pec tively . 

k 

Let V acnote the floi/ of commocioir k berween its 
source-sink pair and x^ denote the flow value of commodity 
k on ax*c (i, j) . 

ObiecTive funct ion ; The objective is to maximize the sum 
of flows of bhe uwo commodities and it may be expressed as; 


Ma^:. (v^ + v^) 


(3.1.1) 
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Constraints s I'liere are three diflerenc sets of constraints 
for the problem. 


ri ow ccnserva Gion constraints ; 1‘hese constraints essentially 
represent the fa,cc that flow of each coinuiodits/ is conserved 


at a.ll noaes, encepx source and sink. 

" k . „ „ k 

V if 1 = s 


n 


X 


3=1 




n__ 

3=1 


■31 


{ 

= -V" 


k 


II 1 = 'G 


1 = 1,2s 


.,n (3.1.2) 


'\ 0 otherwise, k = 1,2 


Capacity cons traints s Sum of flows of both the commodi- 
ties over an arc should be less than the capacity of that 
arc, which is expressed as; 

Jk- b, .S ■V'(i?3) e A (3.1.3) 

Ar c-dis noint c ons traints ; Each arc allows flow of at most 
one commodity over it, which may be accounted as; 

^ij ■ "13 = ¥(1,3) £ A (3.1.4) 


Fon-negativicy re strictions ; Lastly there are restric- 
tions that arc f.low values of each commodity should be non- 
negative ; 

> 0, ¥(1,0) £ A, k = 1, 2 (3.1.5) 

fhe arc-dis joint constraints in tins formulation are 
non-linear cons tnaints and usually are more aifficult to 
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hazLdle than the linear constraints, lliis foimulaiion may 
be reduced to a linear integer programr'iing formulai^ion as 

follows ; 


I nteger progr aiiuu ng formulation ; Lee us iniroouce the 0-1 
k 

■variables y and a large positive nui-foer M, Then the 


arc-dis joint cons'craints may be written as: 

^ij “ "’ij ^ ’ ■\/’(ii,j) c A (3.1.6) 

- M y.^ ^0 5 (i, j) e A (3.1.7) 

^IJ ^IJ sA (3.1.8) 

y^^ = 0, 1 , '\/-(i,j) e A, k = 1,2 (3.1.9) 


On arc (i,j) flow of commodity 1 can be positive only 

if y^^ = 1> xvhich follows from (3.1.6), Prom (3-1.7) it 

follows that on the same arc flow of comodi ty 2 caji be 

2 

posi-cive only if y =1. Prom (3.1.8) it follows tha'c 

1 tl 

only one commodity can flow arc (i^j), 

I’or a ne o-;rork having n nodes and m arcs, this formula- 
tion will reqiure ^m variables and 2n-i-4m constraints. It 
may be solved by any one of 'the integer programming 
algorithms, ho-E^ever solution of a moderarely sized problem 
by such an appro onh may be prohibitive because of large 
computer memory requirement. Moreover, cost of obtaining 
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'Che optimal solution may be high because of slow conver- 
gence properties of these algorithms. 

Hence there , is^need for methods iliich overcome these 

A 

liLaitabions , ihc I ranch and bound algorithmj specially 
written for this class of problems, ma^ be applied to it 
successfully, ne bails of the branch and bound algorithm 
are given beloir. 

3.2 BRANCH Mu BOUND ABGORITHM 

I'he branch and bound algorithm seeks an optimal 
assignment of arcs to the two commodities. In ilns process 
of seeking it implicitly enumerates all possible assignments. 
It performs it systematically and aims to do it with least 
computations. 

In the beginning all arcs allow boch the commodities 
to flow over it simultaneously. An appropriate arc is 
chosen as the ba.sis of branching ana two subproblems are 
generated. One subproblem has the chosen axe assigned to 
commodity 1 and second to commodity 2. Reasibility of the 
subproblems is checked and subproblems found infeasible 
are added to 'blie list of infeasible sub problems. A 
subproblem is 'caken out of this list and more subproblems 
are generated in similar fashion. Pruning and fathoming 
actively cut down the ntunher of active subproblems. The 
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algoritlim bernnnooes when list of infoasihle suoproblems 
IS eiapty, llie strategies which chars cf era zc the Draach 
and bouna algorithm are aescribed below: 

U pp e r bo ui iding s Jr a c e gy ; lech sub problem has a uniciue 
subdivision 01 avcs into three following subsets: 

: set of arcs allowing only comciodioy 1 to flo”- over it 

S2 • set of s,rc : allox^ing onljr cormoditj/- 2 to flow over it 

Sq : Set of arcs allowing both the corj^ioclities to flow 

over ic sirnlcaneously o 

To obcan due upper bound, for a subproblem solve 

Single commocLitv i'-tciimal flow problem Iron to over 

uhe netvrork w'^nch coruprises of arcs S-|USq. Repea-t che 

same for second coiuiiocLity over the netTrorh comprisnv of 

arcs S2U&Q* fe- ard v denote bhe maHiial Ilow’' values 

respectively. ihen the value of che u_jpcr bound for che 

1 2 

sub problem is (v + v ) , 

This IS a logical upper bound n c’le sense that any 
fur Cher assignment of unassigned arcs to the comaioaities 
shal 3 - reduce it^ value because manirael Clow problems arc 
solvea ovex reda'ccd networks. 

hr a nching s era cc gy ; Brajichiiig is done at a sub pro ole m. 
Using the arc flow values of the sub problem, winch are bhe 
values obtained \nnle deteriinning xlio upper bound foi' che 
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sulprolleiiis ail o,rc is cliosen which lies si^aultane ous flow 
of both commouiGies and. raasirauxa sm of flox s among such 
arcs. j. ns arc se_ves as a basis of brancljin^ at the sub- 
pioblen and biro new subpro olei'is are genGicteu. One sub- 
problem has bliG chosen arc assigned co commoai'cy 1 and 
second subproblexi has the chosen arc assi,_,ned to commodicy 2. 

fhe 3l__,oiibl-a maintains a lisc of s ub problems. The 
tt/j-Q new subpro.l" is generated at 'chis seep are cbecLed. 

If found infeasible, they are added to the end of the list 
in the order of ircreesing upper bounds, ifeasible subpi’oblems 
are rsed to impj. ovo the incumbent. 

Sea rching stra b ej^r ; This strategy fines our the sub problem 
at which branching is done. Ilie last sub problem of die 
list of sub problems is selected to be blanched aG. fhis 
results in wha'c coj.nron3_y known as 'depuh first search 
strategy ' . 

3-2.1 St epwi se script ion of the algo ri thm 

Pollowiij-g IS the ste Jwise descnpGion of rhe branch 
and bound algoriGjia for two-commouity ec. -.-disjoint iimm>cimal 
flow problem; 
let 

L = / s < be b]ie set of active infea^sible sub problems 

I denote the value of the incumbent 

Ij(s^) denote the upper bound of the subproblem s^. 
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Step 0 


Soep 1 


St ep 2 


Seep 3 


: Set 1=0 and L to be enpfcy. 

Solvu clie single cormodi fcy aa:cir-ia.l flow problem 

for comr.odity ] over networL L^1 TjA)= Repeat 

1 2 

me ea, *3 for second comodic^r. peG v e^iO v 
denooe i-lie maximum flo-^’- val>.es obtained 
respe c tively. 

If ci'c floxr values axe feasible i.e., ao arc 
hevij.1, , fcuc commodities flouixig over it^ stop, 
it J.S tie optimal flo-jj otlie^'vTim leo subproblem 
Sq represent these flow valros tTirh upper bouna 
as ( '7 -1- v^) . 

Add 10 list L. Go 'CO seep 1. 

; Choc'' che list L. If it is empty go to s'tep 6, 
othennse 'cake ouc che last subprobloi-i of che 
list. let it be UpOatr u'-'e list i, and go 

'CO Seep 2. 

; If upjer bound of s^ is less than the incunibeii'c 
then '^riuie this subproblem an'-i. go to step 1, 
othermee using che arc flO’;j values of the sub- 
pro bl.em, decsrmine "che arc v^dncli has flo-ir of both 
comi'iodi tio s over it ana rQaxi'iUm sum of the flows 
among such arcs. Go to step 3. 

: Geuea.ace two subproblcms s^ and Subproblem 

s^j has che arc selected in step 2 assignee to 
coixiioci oA/' 1 and subproblem has the selecced 
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arc acsi^iied to coinmodity 2. ijetei'iux-ie chs upper 
■bound : lor "both su'bproblens and go to step 4. 

Ste; 4 2 Clieclc r.L suLproblei'i lias loasrble ilou values. 

II so ; to step 5 . 

Si'mlcx’ly c'iieck for sub problem. 

xrfear-'rJ.le subproblems to 'die exid of list L ui 

mere a J rag order of the up^per bounc. 

Go CO srep 1. 

Step 5 s If upper boimd of the subproblem rs greater ohan 
the -value of the incumbent, ohen change 'che 
incumbent value, the new value equal to the upper 
boimc of the sub problem. 

Go to step 4 • 

Step 6 ; Inc LU bent value is the maximal floii- vclae. Scop. 

3.2.2 

A Small 11 - u work flow- problem is solved depi-cting steps 

of the branch ano. bound algorithm developed in this section. 



NOTE; The arc number and arc capaxitj ore mentioned over 

each omc as circled ana uncircled numbers respectively. 
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Il TRAllO F 
Step 0 ; 


Seep 1 ; 

Step 2 s 

Step 3 ; 


Si^ep 4 ; 


I lSRATION 

S^ep 1 ; 


1 

Sec 1=0 aJid I to be emp'cy. Solutions of tlie tuo 
S.c. iienimal flow problems yields = 90 and. 

Vg = 75. Plo\ values are infeasible, let these 
flof- values correspona to subpi’oblen (l)''jrith 
upper bound as (90+75) = 165. Hxls sub problem is 
added to list L. 

Sujproblem (l) is taken our of the lisc L. 

Go to step 2. 

Upperbouiid of the sub problem is greater chan 
incumbent. Prom the flow voAvus at tins srb problem 
the arc selected as o basis of branching is arc 7. 
I'wo uxeu sub problems labelled a.s (2) and (3) are 
generaced. Subproblem (2) has arc 7 a,ssigned to 
coiomoc'ity 1 ana su.bproblem ( ^) has arc 7 assigned 
to cOiiLiouity 2. Upper bounds aie founa as 
U(2) = 155, U(3) = 125. Go GO step 4. 

Both Un subproblems are found co be infeasible 
ana arc moved to list L in ox’oer of increasing 
upper bound. Go co step 1. 

Subpi-oblem with label as (2) is taken out of list 
I. Go to step 2. 
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Step 2 : 

Step 3 : 

Step 4 : 

Il'ERAT IOIT 
Step 1 ; 

St ep 2 ; 

Step 3 ; 

Step 4 s 

Soep 5 ; 

IlERATION 
Seep 1 ; 


Upper ‘bound of tMs subpx'o'bleiii lo greeeer than the 
inciual^ent . Arc 6 xs selected as che "basis of 
■braj-iCnng, &o to step 5. 

Subprobleiii with label (4) an.d (3) are genera tea 
\rith upper bounds as U(4) = l40 , U(5) -- 125 
respecuively. 

Both i^ne s ubprobleios are infeasible and hence 
are ricred to list L in oi'der of increasing upper 
bound, G-o bo step 1, 

-3 

Subproblei (4) is taken out Ii'cn list 1. 

Go to c cep 2 . 

Upper bound of this subprobleu is greater xhan 
incuKibent, Arc 5 is selected as a basis of 
br ail c lung . 

As the result of brandling subproblem (6) and (7) 
are ge' ora ted with upper bounds as U(6) = 120, 

U(7) = 120. G-o to step 4. 

Both une subproblems generated are feasible hence 
are not adaed to L. Go to step 5. 

I = 120, Go to step 1. 

.1 

Subproblem (5) is taken out of list L. Go to 
step 2. 
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step 2 ; Upper bound of the sub problem is greai^er than 
incumbent. Arc 5 is selecced as a basis of 
brancliL/ig. Go to step 3. 

Step 3 I Subproblem (8) and (9) are genera ced as a resull of 
bran Cling ui th upper bounds as U(8) = 120 5 
U(9) = 110. 

Stexe 4 t Both uhc sub problems are feesible. Go to step 5. 

Step 5 i 1 - 120 s Go to step 1. 


Ihis \7a.j sieps are to he folloueu. The branch ana 
bound tree developed by the algorithm for tins problem is 
as follows: 
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5.3 HEURIS'I'IC MLxHOUS 

I'he ■'oraiicli ar.d bouna metlioa developed for tlie xxro- 

comuiodi by arc-uxs ooxnt mazimal flow problem can solve a 

I'loderabely sxzed problem in a reasonable amoixax of Ixme., 

flow 

Solutions of a, large network^ problem by tins approach may 
be prohibitive eve lo large computer meiiory requirement and 
"Ghus a need aric.-s for the methods to solve large problems, 
ilso there may to situations where only a good solution is 
required and search for the optimal solucDon may nou be 
necessary. Moreover, a tradeoff between computer resources, 
i.e., time and moiaory, ana the optimahicy of the soluxion 
may suggest that recourse may be taken bO the heuristic 
methods. 

In this section we shall uevelop t-^/o heurj.sGic methods 
which are intuitively appealing and an-e eeny co implement. 

5.3.1 Heuristi c jL 

fhe laea, behind this hearistic is very simple and yet 

eery effective. two single commodity mamml flow problems 

11 2 ^ 

are solved be tureen s uo t and s to t" rospecxively . 
fLe arc which has flow of both the comnodit les and lip,s 
maximum sum of flows among such arcs is selected, decrease 
in the value of the maximal flow for each commodixy is 
calculaxed if tie selected arc is removed from xhe netx^ork. 
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Pinally, the selocL'ed arc is assigned to che coiinodity for 
which decrease in ihe maximal flow is larger. Plow values 
are changea acccrcangly and the steps ero repeated, When no 
arc IS left wiucli has flow of both co oinio di fci es over it, 
the lie t hod terminates. 

The detailed step by step procedure is given below; 


itotatiois 


; Set of ai’cs assigned to commodicy 1. 

S 2 ; Se b of arcs assigned to commodity 2. 

Sq ; Set of una,ssigned arcs, 

0 ; Pmpty sot 

Step 0 ; Set 5. = S 2 = 9^ 

and ^0 “ • 

Go bo 3tep 1. 


Step 1 ; Solve s.c. max, flo\j' pi-oblem from s^ to over 

G(K,A) , Let X = \x i denobe ohe arc ilow vector 

I ^ 

and V" the maximal flow value. 

lepeob -he same for commodity 2. 

2 ( 2 1 2 

let X = \^iji ^ denote bhe respective flow 
values. Go to steo 2. 


1 

Sbep 2 ; Using the arc flow values X and select the 

In P *1 ^ 

arc (u,l) such bhat X^ 'g 0, Zj^'^0,and 

IS miaximum among all such arcs. If no such arc is 

found go bo sbep 7j otherwise go to step 5* 
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Step 3 


Step 4 

Seep 5 

Step 6 

Step 7 


Solve the s.c, max. flow problem from to fc"' 

over C(]ls,A’); where A' = SqUS^ - (k,!), let 
1 

and "v:^ denote the arc floij- vectoi- and maximal flow 

values respectively. Similarly do it for second 

2 2 

commodity from s to t over Gr(il,A") where 

2 2 

A" = S^hS^ - (kjl). let X and v^ be the respective 
flow V. Ires. Go to step 4. 



w O O 

If (v - v^) ^ (v^ - v^) ,go to saep 3; otherwise 
go to step 6o 

Assign the arc (k,!) to commodity 1. 

= S^ -1- (k,l) 

Sq = Sq - (k,l) 

2 2 

Z = Z,. 

2 2 

V = V, 

TS 

Go to step 2. 

Assign the arc (k,l) to co’Pjjoait 3 ’- 2. 

S2 “ i-'p d- (kjl) 


^0 “ ''"U 

yl -1 

A - 

1 1 

V = v_ 


(k,l) 


Go to a tep 2. 


1 p 1 p 

be floi; values X , Z , v and v are the opximal 


flow values. Stop. 
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3 . 3.2 Heuristic. -2 

‘the unaerlying thumb rule in b'ns heuristic is that 
augmenting the flow along maximal augmenting loath yields more 
flow, definition of the maximal augmenting path ejici an 
algoribhm bo find such a path is given in Appendij>. A. 
Procedure starts \ aebermining maximal jugmenbing pabh 
from s to t and siiiiilarly from s uo . Jlo?r is 
augmented along clie path which has more augmeni:ing capacity. 
All the arcs thob lie on the path selecrec’ for augmentation 
are assigned bo bhe commodity which flows over it. Oonbrary 
to heuristic -1 more than one arcs are assigned to a commodity 
in one icera.tion. Hith the changed flof^ values^ maximal 
augmenting pa/bhs are debermined for each commodity and path 
■\rith more augmenbing capacity is augmented by the respective 

commodity. I'his step is repeated irnibil no augmenting path 

11 2 2 
exists from s uo t and from s to t . 

A detailed step by step procedure 3 s given below: 

The sets S2 a^'^d Sq have the sam.e meaning as in the 

description of heuristic 1. 

Step 0 s Set = S2 = Sq = A. 

Go to step 1 . 

Step 1 ; rete'TiLine iiiaximal augmenting pabh from. s~ to t^ 
over bhe network which comprises of arcs S^USq. 
leb a-i denote the augmenting capacity of this path. 
Go "fco B 2« 
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Step 2 


Step 3 


Step 4 


Step 5 


Sbep 6 


peterimne maximal augmenting patli from to t^ 
over ere network which comprised of am’cs S^USq. 
let denote the augmenting ccpacicy of rliio path. 
Go to oiep 3<. 

If a-^ = 0 and a 2 = 0, go to ouep 6,- otherwise 
detenaine xhe larger of rhe ^wo numbers. 

If a-, a^ go to step 4j Ouhemrise go to vsxep 5. 

Augmeiio the flow on the maieimal augraen tiiig path 
of coi-niodity 1. let denoee ohe sex of arcs 
lying on this path, fhen, 

51 = S3_tJS^ 

^0 ^0“^3 
Go to step 1. 

Augment the flow on the maxiinal augmenuing path of 
commodity 2. let S^ denote ehe set of arcs lying 
on this path. Then 

52 = S2US^ 

^0 ^ 0“^3 

Go to juep 1. 

The flo^r values are optimal flovj values. 

Scop. 
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3 , 4 COMPUl’ API 01 ■ IJj PDRPORJIilTCE 

The brand and bouna aligorithni anc bo Lli die lieurisricSy 
aescrxbed in pro\rioue sections of xlns chapters ■'jere coaed 
in F0R'fRAN"IV for rEC-1090 tine --sharing computer s/stem, 

I'iie compuber orop.rams were tested over o. large numoer of 
random networks generated as folloirs' 

I'he node-node adjacent matrix, iduch specifies a 
nebwork, was gcuierabed randomly, for each position of the 
adjacent matrix a. random number, uniformlj'" distributed 
over interval 0 to 1, was generated. If the number was 
found to be less cliaxi a prespeciiied ii'inbeL, then 1 was 
placed in than position, otherwise 0 was placed there. 

Some other considerations were taken inio an count which 
insured that tho network would be connocced graph and there 
would bo at most one directed arc between any pair of 
nones. PositiO'is for sources and sinks were generated in 
a similan mann.-._ , Par ammeters of a raiicOihLjr generated 
network are number of nodes in the ne tmrk and its density. 
¥ith number of nodes fixed, higher dei^-sicy means more 
arcs. The relationship becween number oi Y>OcteS(n), number 
of arcs (m) an-- density (d) is 
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Ac the end o„ each computatzoiiaTi rua two caiae 'values 
are available - CPU time and elapsed bime. CFJ (Central 
Processing Unit) rime pro'vides the time chat CPU spent to 
complete the compu'cations , vrhereas elapsed time r'movides 
the time curing -jliich the programme \ra.s m che state oP 
execution. JJI]G~10S0 system is a multi^iro gramming system 
ana a programme remains in the stale ol execution for a 
relatively long time ifhile CPU processes it only inxermit 
tantly. fherefo'ce elpased time is generally more than the 
CPU time depentang upon the load of the systeiii. I'his is 
the reason that Cj U time is considered the right estimate 
of computationaJ. time. 

Pwo typec oC eZj_ ei^imen'cs uere concucted. First ty^e to measure 
the computational time of the computer programs of all tbe 
algorithms and second type to measure tbe aeviations of 
the heuristics from the opcimal solution. fhese nere 
conuuccea for sc^^eral combinations of the parameters, "with 
number of nodes varying fi’om 10 to lOu ana aensity varying 
from 0.10 to 0.00. Ihe branch and uouau algorithm could 
solve problems having 50 nodes and 200 arcs, trhereas 
heuristics could solve problems Ixaving 100 nodes and 
2000 arcs. 

fo C 0 iiv..uct first type of experiment 10 randoefLy 
generated problcxiS T.reTe solved for each combination of the 
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parsmetei-s by all algoriciiins and CPU cirae nas obsex'ved. 

I’lae results of tlie experiment are tabula Led in fable 5.4.1. 
'three figures in each, square of the tabDe represeiit the 
average CPU cin’e in seconds xo solve Oxae _^xroblera by branch 
and bound algoritbii, heuristic-1 and heuristic- 2 respec- 
fcivel 3 ’', I'o conduct second type of experiment 1C problems 
were solved for each combination of Ghe j^arameters. Results 
of this experix-nont are tabulated in 'table 5.4.2. I’lns 
table proviaes percent problems t/hose riarcimal flow value 
deviated from optimal flow value b^/ 0/, 0 - 5 ^^ 5-10/ ..., 

20~25y^ or morOj when solved by heuristic-1 and heuristic-?. 

to derive more meaningful results from the tables and 
to get the feel of changes in the computational times j 
several graphs wore plotted. In Pig. 5.4.1 the computationa,! 
time IS plotted vs. density for all che edge-dis joint 
algorithms. tins graph shows the variation in times as the 
network becomes more dense. Pigc. 5.4.2, 5.4.5 and 5.4.4 
demonstrate the relationship between computational time, 
density and no..es for each edge- disjoint algorithm. 

5.4.1 P er f orui^a^g ce _ of the branch and Pound Algorith m 

fhe number of feasible solutions for an ai'c-dis joint flow 
problem are proportional to « Thus as the problem size 
increases, the number of feasible solutioiB increase 
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of 

enormously, fven though, the percen'coge/feasxole solutions 
enumera,fced hj ihe oranch and. bound algorithm deci eases 
with increase in problem size, rhe toucl enumeraced solU" 
tions still increase rapidly, I'liis di’awback makes chis 
uneconomical for large problems. However small or moderace 
sized problems a: e solved very efficienrly . If the number 
of arcs 2 X 1 a neu’[7ork are less than 100 rhen it is generally 
solved in less ulian 5 seconds. So we recommend the use of 
the branch and bound algorithm for smadJ. and noe'erale sized 
problems. 

3.4.2 P erf orm axice s of the Heuristics 

normally computational time and accuracy of the 
solution are the criteria which measure the effectiveness 
of the heuristic programmes. Both the heuriscics in this 
work measure ’trell on these scales. Ihoir efficiencj’’ xn 
solving large problems within a reasonable time recommend 
their application. Problems having IGGb arcs are generally 
solved in less bhan. 5 seconds. 

Both che heuristics provide q_uiiie accurace solutions. 
At least 60 porcenr of the solutions pro'/ided by heuristic-l 
were the optimal ones and in 90 percent of the cases 
deviations from optimal solutions were iichin 5 percent, 

Hone of the problems deviated by more rhan 20 . Tieuristic-2 
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also provioes CLi.ij.-ce solucions, luc .lo'c as ^ooc. 

iieuris cic-1. 

Genera-U}- lieuristic~2 solves a jrohlenin less "ciiae 
tlian lieurisbic- 1. Hovever "bhe ailf ercace in tlie tine is 
not .^reat, I’o siui' arize, we recoranen ^lie use of br8j.T.ch 
and bound al^oritlnii for sioall and modeirte sizec. problems, 
and applicacio.-' of lieurisbic-1 for Isx : o sized xjroblems. 



j-'abla 3.4.1; Coiaoaua aional times for e cl ^e-dis 301111 * algo- 
n tins (I'lmes are in seconds) 


Po . 


jjensity 

o.iw 0.15 

of Nodes 


10 


15 ■ 0.02 

0.02 

0.02 

20 l.lo 0.05 

0.0^ 0.03 

U.03 0,04 

25 C.l4 0.05 

0.04 0.04 

0,04 0.06 

30 „ « 

C.'-^o 0.09 

L . l 7 0.09 

50 

0,25 0.66 

0.53 0.46 

80 - - 

1,67 1.59 

L.o3 1.63 


1.54 2.55 

2.35 3.09 


0.20 

0.30 

0.50 

0.80 

0.01 

U.G2 

0.05 

0.25 

0.01 

0.02 

0.02 

0.07 

O.ol 

C,02 

0.02 

0.04 

0,02 

0 , 05 

0.45 

0.76 

0.02 

0.02 

0.14 

0.40 

0.02 

0.03 

0.07 

0.13 

0.05 

0.12 

0.53 

2.30 

0.04 

i. ( 0 

0.45 

0.79 

0.04 

0.08 

0.15 

C. 22 

O.IC 

5.03 

6.53 


0.09 

0.18 

0.95 

2.53 

0.09 

0.12 

0.24 

0.58 

0.13 

0.57 

1.02 

3.05 

0.14 

1 . 24 

0.47 

0.89 

0.G6 

1.24 

2. 55 

— 

0.6A 

l.v7 

2.31 


2.10 

3.88 

: 

: 

2.53 

4. 18 

— 

““ 

3.45 



““ 

4.55 





100 
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'a"ole 3.4.2; Perf oriao^ico oX edge- disjoint heuristics 


Ileurisiic- 1 


He uric, bic- 2 


To . of 

jjeii 

Percent proMems with devia 

■'o-rceno 

Ox L 

3 0 jloas with 

lodes 

sicy 

tion iroi' 0 p'u'i-^al as; 

oLerici 

on Xrom o_ oimaM a-S; 



0 5 10 15 20 25 More 

0 5 

10 

15 2u 25 More 

lu 

G * 2 

100 

1Cl> 



10 

0 o 'G 

10 0 

10 L 



10 

0.4 

100 

ICO 



10 

0.5 

100 




10 

0.8 

60 30 

90 

10 


15 

0.2 

100 

loO 



15 

0 . 3 

lOt 

So 


10 

15 

0. ^ 

80 10 10 

SC 10 

10 


15 

0.5 

SO 1C 

100 



15 

0.0 

60 40 

30 50 

10 

10 

20 

0.1 

100 

ICG 



20 

0.2 

90 lo 

SO 


10 

20 

0.3 

90 10 

100 



20 

0.4 

100 

70 

10 

10 10 

20 

0.5 

80 20 

80 10 


10 

25 

0.1 

100 

ICO 



25 

0.2 

90 10 

CO 10 

10 


25 

0.3 

70 20 10 

40 y C 

10 

10 10 

25 

0./1 

80 20 

8c 10 


10 

30 

0.1 

100 

lOo 



30 

0.2 

90 10 

90 Ic 
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CHAPTER lY 


MA2IMli PLOT’ 1/ NOLE-DISJOIHT TyO-GOIEIOLITY HET¥ORiiS 


This cha' consicers another vaTriation of che tij-o- 
cotiiniodity maxims'.', flow prohlem. Additional constrs.ints are 
adclea to the general problem that no node allows simulta- 
neous flow of 'two commodities. The resulxing flow x^at tern 
IS node~dis ooina m ihe sense that both commodities flow 
over disjoint set of nodes. 

A result which is immediately avs-ilable from the study 
of node-ais joinc flow is that it is arc-dis joint also, which 
means that node-dis joint flow is a furi.her restriction of 
arc-disjoint flo-% Moreover all node-dis joint flow problems 
may be converted to arc-dis joint flow problems by a simple 
transformation cescribed below. 

Eor each node which hs.s indegree and outdegree 
strictly grester than one, and one ne-''^ node, which is 
corjjiected wixh ole former by a directed ai’C of large capacity 
such that all in.i.Ciin.ng nodes are incident upon the former 
and all outgoing on the latter. Other nodes need not be 
ais curbed. Thus the transformed network, in general, will 
have more nodes and more arcs than the original network. 
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Thxs transformr bion is best understooa a simple illiostra' 
tion given belows 



One nay 1 j tempted co apply tlae oecliniques of arc- 


clis joint flow lino lien to the transformed network. However, 
a little reflect:' on will reveal that ib may not be (desirable 
to cLo so, because ceclmiques speciallj'’ developed fox' noae- 
dis joint flow in-’l solve a problem in xiucii less time chan 
the applicaoion of arc-dis joint flow ceclmiques to the 
transformed network. Reason is simple - number of nodes afe 
geiierallj^ less chan 'che number of arcs anci an algorithm 
whose basis is assignment of nodes will yield the solucion 
in much less 'cire than the algorithm u'hobe basis is assign“ 
ment of arcs. 

Ihis chapter follows the same pstv^ern as thok of 
Chapter 3 . Section 1 discusses the rua t jicma'Cical programming 
formulation of L.ue hro-commodi'ty node-dis joint maximal flow 
problem and an mrceger programming formulation for the same. 

A branch and boux u algorithm developed for this problem is 
discussed in ,3eccj.on 2. Section 3 describes the txfo heuris- 
tic approaches lor the same problem and section 4 presents 
the computationcT performances of the algorithms developed 
in previous sections. 
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4.1 MAx'HSMllI CiiL P-iOGRAJ'S^IFG PORMLA'-HOil 

Ob'ieclxve Piu ciiioii : The objective is go raaiiimis’e the sum 
01 floirs 01 the Giro commodities which iiiey be expressed asi 

kn: . ) ( 4 . 1 • 1 ) 

C onsiram bs : ‘Ucj e are ihree differeiib sets of constraints 

for this pro bl eii . 


Plow cons erva'Gi Oil co nstraints ; These cons'crainis represent 
ihe fact bhac flow of each commodity is conservea at all 

nodes exce pt scarce and sink. This nay be written as 

k . k 

V II 1 = s 


n 


3=1 


k 

"13 


X 


n 


31 


3=1 


=^-v^ if 1 = , 1=1, 2s . . . k - 1,2 


i,jD otherwise 


( 4 . 1 . 2 ) 


Capacity const rsinbs 5 Sum of flows cf both the commodities 
oyer an arc should be less than the capcnity of that arc; 


X e A (4.1.3) 

lj=l ' iJ 


K ode-dis .ioint_ c onr. trainbs ; Each node allows at most one 
commoaity to flOG through it. It may be cccounced as; 

ii- 2 

!•_ , * aT- 1-., =0 , 1 = 1, 2, n (4.1.4) 

3=1 3-1 ^ 

E on --ne gat lyit y r e s t ri c t ion 3 ; Lastly .here are restrictions 
blast arc flows of each commodity should be non-negabive ; 

0, '^(1,3) e A, k = 1,2 


( 4 . 1 . 5 ) 
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The noae-ois Qoini; constraints of ohis formula'cion are 
non-linear in nature and difficult to handle. They may he 
transformed to linear constraints as follows ^ 

]£ 

I nteger pro^ ra i'Hii ift formulation ; let che 0-1 variables y^ 
and a large positive number M be introduced. Then the node- 
dis^oint cons'craiii cs may be writuen as: 


X ~ yi 4 0, i = 1,2, . 
3=1 “ ^ 

S - M y ^ 0, i = 1, 2, . 

3=1 ^ 

yj^ + y^ 3 , 1 = 1,2,. 

y^ = 0,3-, 1 = 1,... 


« ;n 

(4.1.6) 

. jll 

(4.1.7) 


(4.1.8) 

n j k = 1 , 2 

(4.1.9) 


Prom conscrainbs in (4.1.6) it follows tha'c flow of 
comifloaity 1 through any none Ni can he positive only if 
y^ = 1. Constraints in (4.1.7) impose the same restriction 
for commodity 2. I'rom constraints in (<.1.8) it follows 
tha,t at iH0st oxie commodity can flow ilirou-gh node hi. 

If a network which has n nodes ano m arcs and is 
formulated in ^ i.e xfay described above, chon it will require 
(2n+2m) variables and ( 5n+m) constrainas. *This formulation 
can be solvea 63- any one of the integer program iing algorithms 
however solution of a large problem by chis approach may be 
prohibitive due jo large number of variables and constraints. 
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Moreover cost ox otoaming tlie optimal solution iia,y be high 
because of slor convergence properties of these algorithms. 

i’he branch and bouna methoo, specially ^rntlen for 
this class of problems overcomes these lii^n tat ions and solves 
a moae'vtely l&rje problem in a reasons de amount of time. 

4.2 BRAHCH iiJTu tOUhn AlGORI'lHM 

The branch caici bound algorithm seeks an optiria.1 sub- 
division of iiooes to the two commodities. Initially all 
nodes allow simuluaneous flow of both the commodities through 
it. An appropriate node is selected and tuo sub problems 
are generated, one subproblem has the selected none assigned 
to commodiTy 1 ana the second subproblem has olie selected 
none assigned to commodity 2. Feasibility of the subproblems 
is checked and infeasible ones are added to a list of 
subproblems. Ah appropriate subproblem is taken out of this 
list ana more subproblems are generaied. Pruning and fathoming 
actively cut do-n che number of subpioblems . The algorithm 
terminaces w]'eii list of infeasible subproblems is empty. 
Following are the straregies which characterize the brand 
and bound algoriulim. 

Upper bo und.ing s tr ategy ; dach subproblem has a unique sub- 
division of nouec mto three following subsets; 



; Set 01 noaes allowing only coioi’iodifcy 1 to flow 
through it. 

S 2 s Set ol nodes allowing only coauocity 2 to flow- 
through it. 

Sq i Set of iiodes allowing both the coomodities to flow 
through it. 

'i’o obfcaiii gI e upper bound at a sub problem solwe single 

commodity manmu-i floxf problem for conmodity 1 over the 

network comprising of xhe nodes S^^USq and arcs bha'c connect 

these nodes. Let v^ represent the manmu-i flow volue. 

Pi.epeat the same for commodity 2 over the network comprising 

2 

of nodes 321130 al’u arcs connecting these nodes, let v 

represent the maximum floxf value. Then the value of the 

1 2 

upper bound for the sub problem is (v - 1 - v ) . 

Branching strat e gy ; Branching 13 doae an a sub problem. 
Usin_, the arc flow values of the subproblems s ^diich are the 
values obtained while dexermining the '^pper bound lor it, a 
node IS selected irhich allows simultaneous flow of both the 
commodities and .laximum sum of flows among such nodes. This 
node serves as a basis of branching and at a result two 
subproblems are ,ene raxed. One subproblem has the selected 
node assigned bo co.imodity 1 and xhe second subproblem has 
the selected node assigned to commodity 2. 
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The algoi'ZLolim mainxaixis a list ox sub problems . The 
two new subprobl^uTS generated at this step 8,re checked. If 
found infeasiblSj ohejr axe addea to the list in the order 
of increasing upcer bound, feasible sabproblems are used 
to improve the incumbent, 

Ss orching stra x e.cy ; This straxegjr finds out the subproolem 
a,t which branching is done. The last sub problem of xhe 
list of problei'is is selected for further branching. 

4*2.1 Stepwise _de script ion of the algo r ithm 

Following IS the stepicise description of the branch 
and bound algorithm for two-commodity node-dis joint maximal 
flow problem. 

ITota-tions are same as described in section 3.2.1. 

Step 0 ; Set 1=0 and 1 to be empty. 

Solve the single commodity mapcimal flow problem 
for commodity 1 over the netirork Qr{l^,j±) ^ Eepeat 
the sane for commodity 2. Let v^ and v^ respec- 
tively denote the maximum flo-^.r values obtained. 

If arc flow values are feasible, i.e., no node is 
having two commodities flowing through it 
simul caiieously, then stop as it is the optimal 
flowj Ouherwise let subproblem Sq represent these 
flow values with upper bound as (v + v ) . 

Add Sq to L and go to step 1. 
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Step 1 ; 


Step 2 : 


Seep 3 % 


Step 4 s 


Step 5 ; 


Step 6 ; 


Check tke lisx L. If it is enipty go to seep 6? 
otlierTTise bake oux ehe las'c suoprobleu of the list, 
let lu he s^. Update the list 1 and go co step 2. 

If up'er hound of is less than the incumbent 
'cneii 03’une this sub problem anci go to step 1; 
oxhemise using the flou va,lues of the subproblem 
detersiine the node which has ilo^r of both the 
commodi cies over ic and manimuiii sum oi ilof^s among 
such nodes. Gro to seep 3. 

Generate two sub problems s^^ ana Subproblem 

has xhe node selected in svep 2 assigned to 
commodity 1 and subproblem b has che selected 
noao assigned to commoaicy 2. jet ermine the upper 

bouids for both subprcblems ana go to step 4 . 

Check if sub problem has feasible flow values. 

If so j go to step 5 . 

SiiTilarly check for subproblcm ^.dd the 

infeasible subproblems to Che end of list L in 
increasing oxaer of the upper bounds. Go to step 1. 
If upper bouna of the subproblem is greater than 
the value of the incumbent, then chaiige the 
incujabent value, the new value equal to the upper 
bouno of the subproblem. Go to step 4. 

Incumbent value is the maxi.ial flow value. 

Stop. 
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4.3 HBURIMIC ISl'IIGuS 

ihe braxLch and bound method developed for the two 
commodity node~di& joint maximal flow problem can solve 
proDlems having 50 nodes and 300 arcs. Ilotrever, for larger 
problem size exponentially increasing computational bime and 
large memory I’ecinrement prohibib the application of tins 
approach. A recourse is sought throu^i heuristics to solve 
such problems, two heuristics vjhich are along the same lines 
as their edge-dis joint counterparts axe proposed. Iheir 
ae balled discuscioxi follows now. 

4,3.1 H euristi c"! 

fwo single commodity maximal floi7 problems are solved 
from s to 1 and s " to t respeccively . fhe node if/hich ha.s 
flow of both blie commooities through it and has maximum tobal 
flow among such nodes is chosen, necrease in bhe value of 
the maximal flOTr for each commodicy is calculated if the 
chosen node and arcs incident upon it are removed from the 
network. (The chosen node is finally assigxied to the commodity 
for which decrease in the maximal flo^ij- is larger Hith the 
changed flow values, which occur as a result of die assign- 
ment, another noc.o is selected employing the same criteria. 

Ihe method torminates when no node is left which has flow of 
both commodibies through it. 

Pollowing step by step procedure my be suggested for 


this methodology. 
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HOI'ATIOjS 

: Set of roaes assignea to coimj.ioc.it 3’- 1, 

Sg i Set of iiJ^des assigned to coiimocLi 1:3^ 2. 

Sq ; Set of unossigned nodes. 


Step 0 


Step 1 


Seep 2 


Step 3 


Sc b Sj = = P, and Sq = A. 

Gro to Sv^ep 1. 

Solve s.c. max, flow problem fro-ii s~ to t^ over 
G-(FjPi). Let denote tlie arc flow vector 

and v^' tlie maximal flow value. 

2 f 2 ' 

Repeat tlie same for second commodity. Let }[' = 

f J 

2 

denouG the arc flow vector sjiC v the rnaxrnial flow 
value . 

G-o to step 2. 

-1 p 

Using ulie arc flow veccor and Z select the node 

-- T ^ - p 

n such that ^ x^ , 0, x .. X 0 and 
n , n 

( X ■ + ^ X ) IS maximum among such nodes. 

3=1 3=1 

If no such node is found go to step Tj otherwise 
go to seep 5 , 

-j U 

Solve uhe s. c, max, flow prohlem from s~ to t"^ over 

the no Gwork which comprises of nodes S-^USq - WqI 

and arcs connecting these nodes. Let 2^ and v^ 

denote the arc flow vector and maximal flow value 

respectively. Similarly do it for second commoaity 
2 2 

s "CO t over the network which comprises of 


from 
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nodes S^USq “ and ancs connecting these xiodes. 

2 2 

Let Olid V denote the respective flow values, 
s s ^ 

Go to etep A . 

Step 4 ; If (v^ - v^) JX (v^ - v^) go oo step 5, 

52 5S1 

otherw-i .e go to step 6. 


Step 5 : 


Assign the node n^ to coimioca.t5^ 1, 


Snap 7 


^1 

^0 

.2 


= V. 


^1 po'^ 

‘^0 “ r'^o) 

..2 

,2 


35 


S'oep 6 j 


Go to step 2. 


Assign the node n to comaiodiiiy 2. 


^2 - 
So = 
x' = 


1 


V 


= V. 


^2 'Po^ 
Sq 
-G- 
1 


r ^ 


Go to step 2. 

j_ 2 1 2 

The flow values X , v and v are feasible; 

hence stop. 


4.3.2 Heurist ic 

Proceaure starts by determining uazirrial augmenting 

11 2 2 
path from s to a and similarly from s to t . Flow is 

augmented along the paih which has more augmenting capacity. 

All the nodes tba t lie on the path selecrea for augmentation 
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are assigned, bo 'che comnodity -wiiicli flotrs over ic, Witli 

the changed flotr values this basic step is repeatea again. 

1/hen no augmencing path with positive augmenting capacity 
11 2 2 

IS left from s vo t and from s to t , che proceouro 
terminates. 

The sets S 2 and Sq have the same meaning as in the 
description of lieunstic-l in section 4.2,1. 

Step 0 ; Set S-, = = jd, and Sq = A 

Go bo step 1. 

1 1 

Step 1 . determine maicimal augmenting pa^cli from s"^ to t 
over che network which comprises of nodes 
and arcs which connect these nodes, let aj^ 
deno cfc the, augmenting capacity ox this path. 

Go "CO step 2. 

2 I 2 

Step 2 ; he ter XLxie maximal augmenting path from s to t 
over the network which comprises of nodes S 2 U&Q 
aixd arcs connecting these noces. let a 2 denote 
the augmenting capacity of this path. Go to step 5., 

Step 3 ; If a-| = 0 and 12 = 0 go to step 6. Otherwise 

determiJie the larger of the two numbers a^^ and a 2 . 

If s- 2 s go bo step 4? othecrwise go to step 5. 

Step 4 ; Augment uhe flow on the maximal augmenting path 
of commodity 1. let S^ denoce che set of nodes 
l37'ing on the path. Then 
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Step 5 ; 


Step 6 : 


— O-jvJ ^ 

“ "0 “ ^3 
Go CO step 1. 

Aujienc the flox7 on the msniriel aupmencmg path 
of commodity 2. Let deiioce ‘^he set of noaes 
lying of che pexh. 'J?hen 


S 2 = S 

&0 = s 

G'O "G 0 
Prnn c 


G ~ ^3 
seep 1. 

the flow -values 


QJICl S u 0 0 * 


4.4 COMPUl'A'ilOlTiiX PERPORIIANCIi 

I'he branch and bound algorixtuu and both 'che heuristics 
described in previous sections of this chapter 'v’-ere coded 
in EORTRAR-IV for j^LC-1090 time-sharing conputer system. 
These programs >-ere 'tested over a large number of randomly 
generated netwoD’ts. Computational times were measured in 
■terms of CPU iiins. The method to generato random networks 
and justification of CPU time as a measme of tirie perfor- 
mance of algo ri thins is described in sec cion 5.<. 

Two types of experiments ^^ere conoucted. Of first 
type CO measure che computational time of the algorithms 
for several combinations of the parameters and second to 
measure the deviations of the heuristics from the optimal 
solution. The number of nodes varied from 10 to 100 and 
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density varied Trom 0.10 to 0.80. I'lie Irancii one tound 
algorithm can solve problems having 50 nodes ana JOO arcs 
and heuristics can solve problems having 100 nodes and 
2000 arcs. 

To coiiduco firsc ‘cype of esperiiaenc 10 problems uere 
solved fox each combination of blie paname c srs by all 
algorirhms and CHJ time uas observed. rhe results of the 
ezperiuenx are babulac:ed in Table 4.4.1. Three figures 
in each square of blie table represent uhe average CftT time 
in seconds to solve one problem by branch and bound 
algprithm, heuriotic-l ano heuristic-2 respectively. To 
conauct second ej’-pe of experiment 10 problems were solved 
for each combination of the parameters and results are 
tabulated in Table 4.4.2. This table gives pei''centage of 
problems for different noae-density combinations xjhose 
maximal flow vadLue deviated from the o^'timal flow value by 
0 /^, 0-5//, 5-10/, 20 -25 /o or more, when solved b3’- 

heurisxic-1 and heurislic-2. 

Several graphs were plotted to have bectei* under 
standing of the computa <:ional times, figure 4-. 4.1 
demonstrates fie relative pex'formance of no de-dis joint 
algorithms on cinie scale and Pig. 4.4.2, 4.4.3 and 4.4.4 
show the variations in computational ximes as nodes and 
ciensity change. : 

595l9 
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4.4.1 Ba rf q m e loce of Branch and Boun d A l,RoritIiia 

Total nuifo. r of feasible solubions is proj^or bional 
to 111 and soHv. _ over of d. Even tliou^ii the percentage 
of sutprohleuis pejiei’ated decreases wif'i increase in n or 
d, the nu'ifoer o " subyrohlems generated mcroanes rabidly 
^’luh increase jii n or d. However , small sized problems 
are solved veiy efficiently. If number of arcs iii a. 
network are less ohan ICO then there is no significant 
aifference in conpu tational time of solving the x)roblem 
by branch and hour a algorithm or anj" of rhe heui’iscics. 

So application of branch ana bound algorithm is recommenaed 
when problem si^e is small. 

J?roblems ^nth number of arcs less blian 500 are 
generally solvec in less thaii 6 seconds. Hith a further 
increase in nmabor of arcs the computational time 
increases rapiv.13 . It was observed that ''oj doubling the 
number of arcs c f a problem with 500 srcs the computations,! 
time increased oTjc-ntyfold. Still it ic not the computa- 
cional time Tanch restricts the solucion of large sized 
problems, it is the large memory requi 'cement which put a 
limit on tne size of the problem. 

4.4.2 Perform ance o f H eu ristic s 

iwo criteria are generally used to measure the effi-- 
ciencj?' of a heu'cistic method. They ar’e computational 
time and accuranj^ of the solution. 
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Heuns ciC"l solves a large sized ]:rololeiii in quite a 

reasonable tixio. Problems having 100 nooes and 200C arcs 

are generally solved in less than 1 mincLoe. Ic wa.s observed 

quite 

that compuiL'a tiional oime increased^lineaa'’ly i-n-'ch increase 
in problem size. 'IhjLs raetliod also provides quite accurate 
solutioixS. blien the netxTOrk is less dense (d ^0.30), 
tlicri at least SO^ of the problems a.re solved optimally, 
uh fch liicreaseo values of deiisities more problem deviatea 
from opuimal. Jven for large problem sizes at least 
50^ of the pro.lems were solved optimahl-y amd about 80 
of the problems '^rere solved ■^rith &€ mos c 1(^ aevieecion of 
che maximal value provided by the heuriscic from rhat of 
optimal value . 

Heuristic 2 talces much less ‘cine than tliac of heu 
ristic~l but at ohe cost of less accuracy. Hroblems having 
1(j0 nodes and 2000 arcs are generally solvea in less chan 
10 seconds; whicJi is one sixth of the oime co solve 
problem of same size bp heuristic-l, this heunsric also 
exhibits linear trend in computationa]- cime with increasing 
density values. However, the solu cions are not very 
accurate. Accuracy goes on diminislung with increa-se in 
density values. A trade off betTreen comouxational cost and 
accuracy may ueuermine the selection of appropriate 


heuristic. 
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0.01 
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0.04 
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0.07 
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0.09 
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0.13 

0.45 
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0.11 

0.17 

0.26 

0.46 
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6.36 
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0.44 

0 , 60 

1,20 

1 . 64 
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1.73 
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- 
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- 
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le 4.4.2: Pertoriaance of no cle-dxb joint heuristics. 
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CHAPl’ER V 


llAHmOj ArtC-JJISJCIET AITu ITOLE-lilSJOIHl' LrJLlTCOMIOElPY Eioy 

PROBLEMS 

'i'lie pTObler& considered m this chepter are 'ciie exten- 
sions of the probloiiis considered in previous chapters. An 
attempt is made oO develop methods similar to that of t-vj-o- 
commodity flow problems foi the multicoimaoditj?- flow problems. 
Heuristics are exteiiued for multi commcli ay networks, however 
extension of branch and bound algorithms was not attempted 
because branch ana bound algorithms for these problems would 
be computationally very slow. 

In this chapter two heuristics are developed for edge- 
ais joint multi CO viodity maximaJ. flow problem and txTO for node- 
dis joint multi coiiiuodity maximal flow problem. Computational 
performances of che programs, built for these methods, is also 
presented. 

5.1 MAXIMAE ARC-uISJOIHT MOLi'I COffiO rllT PLOW PROBLiM 

I'he problori uO find maximum sum of flows of coEumodities 
in a multi comiQOc-i uy nebworJ:, such that no ai-c is having more 
than one coinmoaity over it, is considered in this secxion. 

fhe problem may be formulated along the sa.me lines as 
the formulation for cwo-commodity arc-dis joint flow problems 
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modifications ttHI de required in the constraints as there 
are q commodities instead of two commodities. The arc- 
dis^joint constraii'XG shall he the most affected one since 
prouuct of arc floiru is zero for each comhination of rhe 
conaiioaities . An iixbeger programming formulation may be 
offered along tlie siurLlah lines as thab of two -commodity 
case . 

Extension of the branch ana bound algorithm for two- 
commoaiiy arc ’■ais;]oint; flow problem to the multicommodity arc- 
disjoint flow problem poses some problems. I'he upper bounding 
criteria will no more remain an efficient criteria because 
all unassigned arcs are assigned to each commodity and as a 
result the upper bound may be as bad as q times any feasible 
solucion. I’hus with increase in number of commodities , the 
upper hounding critei’ia keeps on worsening. This ul’cimately 
leads to large n> foer of suhproblem enumerated and computa- 
tional time IS proportionately large. Moreover^ more flow 
values are stored at each subproblem winch cuts domi the size 
of the problem solved by this approach. 

On accouLic of these reasons no attempt ixas made to 
develop branch a».ad bound algorithm to solve this class of 
problems. 1 ‘he success of the heuristics to solve two commodnty 
flow problem by pi’oviaing solutions very close to the optimal 
and 111 a reasonable amount of computational rime tempted us to 
extend them bo multicommodity cases. The heuristic methoas, 
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wlncji are cescri"o,.d beloWj pioioise a con& icier able saving in 
computat lOiiaH ti ic and proviae veiy gooo soltitions if not 
t iie be st . 

5.1.1 He uristic ^ j. 

q single coi'iaoaicy maximal flov px-oblems are solved; 
each fx'om s'^ to k = 1, 2, .. . q.. ihe arc having Liazini'um 

number of commoai oies flowing over it anc. cribh botal 

flow among such arcs is chosen. Decrease in the value of 
tlie maximal flow for each commodity is calculacea if the 
chosen arc is not included in the network, She chosen arc 
IS fina-ily aSwSigned ro bhe commodity for wliicb aecrease in 
the ma2:imal flow is largest. Plow values are changed 
r? ccordingly and cliis basic step is repea,Ged. ¥hen no arc is 
left which 3,11 ot,'- 3 more than one commodicj' bO flow over it, 
the method terminates. 

I'he Sbep by step procedure of this heuristic follows 

H 0 "v^o 


iT ulAflONS 

s Set of arcs assigned to commodity k. 
Sq ; Set of ujiassigned arcs. 

Seep 0 a Se C k — 1; 2; q.; 

and Sq = A. 

Go to step 1. 
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Step 1 


S'Gep 2 


Step 3 


Step 4 


Step 5 


Solve '^lae single coiimodity maZiKial flow problem 

li! 

from to t^ over G(H,A) for all commodities, 
k = 1,2, let 1 deiiote the arc flow 

L 1 J _j 

vector and v che maximal flow value for coimiodity 
k. Go to step 2. 

1 

Using olie arc flo’.r vectors h = 1,2,...,G, 
select oils arc (a,b) whicb has mammum number of 
commodibies flowing over ib aiio maximum total flow 
aiiiong such arcs. 

If xio sue]' arc is founa go to step 6,* otlierimse go 
to step 3. 

h k 

Solve the s.c. max. flow problem from to t' 
over Q(J1,A^) where = SqUS^ (ab) for k=l, 2 ,... 5 G. 
Let ]C^ ana v^ denote the respective arc flow vector 

srs: JT 

and maximal flow value. Go bo step 4. 
uetennixie the commodity p for ijhich decrease in 
maximal flow value is largeso, i.e., (v^-v^) is 
largest. 

Go bo step 5 

Assign the selected arc 'a, h) bo commodity p. 

Sp = S.^ -! (ab) 

Sq = ■ ( ab) 

= X^, k = 1,2, . .. ,q, kT^p 
v^ = vi^; k = l,2,...,q, k^p 


Go to L uep 2. 
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Step 6 ; The values v"^, k = lj2,...5q. are 

feasible. Stop. 

5.1.3 Heurist ic 2 

The procedrrc starts by determining nasiiiial augmenting 
k k 

paths from s bo b , k = l,...,q. Plow is augmented along the 

pa.ch vdnch has ma-'imal augmenting capacity. All che arcs that 

lie on the path selected for augmentguion are assigned to the 

commodity tj-hicli ''lo’t^s over it. Contrary to heurisuic 1 more 

than one arcs are assigned in one iteration. Hith che changed 

flow Values this basic step is repeated a^gain until no 
p a 0 h T_ V 

augmenting/eiist.^ from s to t , k = 1, 2? . . . ; q. 

Polio-wing IS ohe detailed step by step procedure of 
heuristic-2. 

HOl'P? Sets Sj^ and S^ have the same meaning as in xhe descrip- 
bion of heuristic-1 in 5.1.1. 

Ste p 0 : Set Sj_^ - 0 , k = 1 , 2, . . . , q 
So = A. 

k k 

Sbep 1 : net ermine maximal augmenting pack from s to t 
over ino network which comprises of arcs Sj^USo 
for all k = l, 25 ,... 5 q. let a,^ denote che augmenting 
capac 3 -cy of these paths respeccively . Go to step 2. 
Step 2 ; If av = 0 for k = l, 25 ...,q. Go to step 5, obherwise 
go to seep 3. 
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Step 3 


Seep 4 


Soep 5 


Determine the co iJino di fey p ^u'hici lir„s au^meiicing path 
■^ath largest augmenting capacity^ i.e., a > a, , 
k = 1, 2, . . , , q. . Go CO step 4 • 

Augn^ent che flow on the maximal augmexixing pa'c’a el 

commoui fey j:. Let denote the sot of arcs lying on 

this pacr, I'lien 

S^ = & J j,, 

1 _ 

So = So S^ 

Go to seep 1. 

Print 0 . 1 ? floT/r values and sfeop. 


5.2 MASMAL IIOB>-DISJOIH'r mTiCOMODl'Ii PLOU PP0B13M 

file pro hie 'T feo find, maximum sum of floiirs of commiodi :ies 
in a mulfeicoinmodi by network, such thac .lO node is having more 
fehan one commooicj flowing over it, is considered in cliis 
section. 

Ilafehematical formulafeion of this problem will be along 
similar lines as fox ibs special case of bvo-commodi oy networks. 
Along the similar lines the problem may be formulateo. as 
integer lineal programming problem. 

Ihe reasons, gi'^/en in Sec cion 1, for not developing 
the branch and bcunu algorithm foc'* arc-casjoinx flo-’- problem 
also bold for no L developing xhe same for no de--dis joint 
flow problem. Incfeead two heuristics are proposed. Perfor- 
mance of chese Leuriscics was very encouraging for two- 
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comnoaicy problems and ic is hoped that they xfill al 

do "vrell for nnilbiLCOioiuodixy flow probleius. 


5.2.1 lIeurisbiC‘1 


Sinyle commo^-'icy maJtimal flow problems are solve a from 


s to for k 1,2. 


,q. The node ir’-'ici: al3-ons iiianii/ium 


nuraber of co-miocLi exes bo flo^r thi'ou^^h it and has mariiaum 
total flow a,monp ' uch arcs is selected. ^^ocrease in the 
valuG of the maximal flow foi each con odit_y is csiculateci 
if the selected node and arcs incident upon it are removed, 
from the networm, llie chosen node is finally assi 2 ;ned to 
the commodity .or which decrease in the iiarimal floe value is 
largest. Uich the changed flow values this basic sxep is 
repeated. Ifhen i/’’’ere is no xiode winch alj.ows more than one 
commodity to iIot- clirough i c, the method Terminates, 


ihe S'cep bp step procedure lor tins Eiethod is given 


below; 


N OlAfl OliS 

? Seo of noi-ei assigneu to comiiioait3r k. 
Gq ^ Set of unabGigned nodes. 

Seep 0 a Sec S.]^ = 0 ^ k = i,2,.>.,l; 

J. - 

ant S(_^ = A« 


Go to step 1. 
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Step 1 


Step 2 


Step 3 


Step 4 


Ste p 5 


Step 6 


; Solve uIt^ s.c. max. floxf problem froui 5*^00 t'*'^ over 

1 , 

G-(NfA) k = l^Sy^ocsq,. Lee 1"^ = x'"' aeno te xlae 

I k 

arc flov vector oad v the inarrmaJ. ilo" -^alue o! 
coiimot'i jp k. G-o ro step 2. 

s Dsin^, Xiic rrc flotj vectors 2^, k = 1^2 , select 
the noe'e 11 ^ vhich has ma-ciinmii a umber or coiuaocii ties 
flomirg tlu’ough it and has maximum tov-l Xlou among 
such nodes. If no such iiode ic iou-io.j gc co step 6, 
othen;n.se go to step 3. 

: Solve fTe s.c. mc.x. flow problem Ihom s to t ' over 

the ne t'rork comprising of nodes Sj^USq - n^ and 
arcs conj-'ccting i.-us set of noaec. no it for k = lj2s 
...sQ,. JiC L and v;^^ denote the r>especGive ore flow 
vecoor cod maximal flow value. Go to step 
: heteriii-io the commodity p fo: fT’iT_c]i decrease in the 

maxim 1 Iloir value is larges u , j..e.^ (v^“V^) is 
larges t , 

Go to step 5 o 


AssitX selected node n^ to cominouity p, 

= S.. + n^ 


P 


So = Su 


p 

'0 

,k 


n. 


J 


h = 1, 2 , ... , 1 , k ^ p 


v^ = v^" 5 k = 1, 2j . . . s, 0.5 k ^ P 


Go CO Slop 2. 

k k 

Print flow values X ^ v ^ 


— 3Z1 0. 


s cop 
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5.2.2 He urxs'cic __2 

The procechirc starts by deterffliiixi^ iiasimal augmenting 
Ic ! 

patlnsfro" s to t k = lj2,...s,q.. Flou is augmented along 
the path which has naximun augmenting capacity, iill ahe nodes 
chat lie on the pakli selected, for augment a ion are assigned 
to the commouitjr d ch flows over it. Hence many nodes are 
assigned in one lueration. \/ith the changed How volueSj this 
step lo repeated. Hhen xio augmenting path exiscs from s^^ to 
t ^^5 k = l,2s.c.,Oj uhe procedure terminates. 

The s'cepwioe proceaure now folloxrs. 

NOI'A HON 

Sj^ = Set of noaes assigned bo comriiodit 3 ’- k. 

Sy ; Ser of uiias signed noaes. 

Stop 0 : Set k = 1,2, — ,Q, 

and Sq -- A. Go to step 1. 

k ]c 

Step 1 : net ermine maximal augmenL,ing pauli from e "co x 
over die network wluch comprises of nodes SiHSq 
axia a,rcs connecting this set of nodes, no it for 

kHjlj ,q. Let ai^ denote che augmenoing capacity 

of mavj- lal augmenting path of comnodits'’ k. 

Go to seep 2. 

Step 2 : If = o for all k = l,2,...,q, Go to step 5; 
otlienm-se go to step 3. 
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computational time, vhich was in terms of CPU time, was 
oUserved. As a ror.ult tables were generated winch provide 
average CPU time co solve one problem for aifferent node- 
density combinatioo s. 

ijacli nodO'dcisity combinaoion corresponds to one square 
in the tafole. Ijach square coniains 'chree figures which are 
the computationcJ. cimes to solve proble..i& having 2, 5 and 8 
commoaities respectively, tables 5.3*1 a/'a 5.5.2 present 
the compoutatioual times in seconds for eago-dis joint heuriscic 1 
and heuriSL.ic-2 resjjectively. Computational times of node- 
disjoinr heuristic -1 a,nd heuristic-! are tabulated in 
‘fables 5.3-5 and . .3.4. 

Several graphs were plotted to factlitate better 
understanding of tne computational times. In figure 5.5.1 
computational uiiaes of all the heuristics are plotted against 
densi cj'-. Figure '^.^.2 depicts the rela tioxislnp betxreen 
computational cine end aensity for various number of commo- 
dioies, whereas Pi.a^’® 5.2.3 shows the variation in computa 
ta onal time as density and number of noaos are changed. These 
two graphs are plC'vtea from the table of edge-dis joint 
heuristic-!,- however, saiio xoattem is obseri/ed for other 
heuristics also. 

jjue CO unavailability of exact aJ-gori chms for multi- 
coimnoditj/' disjoint flows, it was no'c possible to measure the 
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de-viai^ions ol o'-., solutions providea clie heuristics 
from chat of o^ ^irs,! ones. Houeverj z’eletive erf or mane os 
of the heuris cico ^ras judged. 

5.3.1 Results ij iC x'lscussip nc 

One obsej’vr oion uiiioh is evident foon Rig. 5.b.l is 
tiaat heuristic 1 takes conciaerably more compucs'cional 
time than heujiGtic -25 he it a arc-di'; jOint floir or node- 
dis^oint floT,*. Iloreover this difference in computational 
times keeps on increasing as problem sise increases. 'This 
drawback of he.’T’istic ~1 is coun cerbal-'iiced by che facx that 
generally it providos better solution,.. The question, that 
which of the t'U is better, does not itc^’-e a sxraiglitforwara 
answer. A troceofX betn-een accuracy' and cost might 
decermine the ^election of appropi’icce heumstic. 

It vas mentioned in chapter A thac rransformacion of 
nod e-dis j oint floxr problems to arc-dis joint flov 
may not be e cOi''Oiiu.cal. Tliis is evidsxit from Rig. 5.3.1 
that node-dis joixic heuristics takes coneiderablj less time 
tlian cheir ar c~dis joino counterparts. Hence, \<[e have been 
right in our decisioii to develop sex^arate algorithms for 
them. 
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ciOiial performance of edge-dis joint 
ic-1. 

j^GlT sity 

o.lO 

FCo of xiodes 

0.15 

0.20 

U.30 

0.50 

0.80 

10 


0.02 

0.02 

0.02 

0.10 


-- 

0.02 

0.03 

0.07 

0.19 

‘ 


0.03 

0.03 

0.09 

0.35 

20 i'.03 

0.03 

0.04 

0.06 

0.14 

0,95 

L .03 

0.04 

0.06 

0,26 

1,64 

9.44 

0.04 

0.05 

0.12 

0.53 

3.72 

11.72 

40 0.12 

0.13 

0.46 

0.75 

1.75 

7.01 

0.23 

1.38 

4.46 

24 . 20 

76.00 

226.83 

0.24 

1,38 

5.44 

37.48 

131.77 

307.62 

60 0,35 

1. 37 

1.48 

2.49 

5.33 

. - 

2.43 

14.93 

29.45 

197.35 

438.17 


7.13 

22.96 

56.10 

314 . 26 

““ 


80 0.90 

2.02 

1. 63 


_ 


30139 

90.99 

257,11 

-- 



y ^ . 7 9 

103.25 

353.32 

* 


•tv*. 

100 :.94 

2.04 





5^i .09 

94.85 

- 

- 

- 

- 

9o.i2 

213.55 
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I'able 5 . 3 . 2 ; 


ueiisxty 
Fo. of Nodes 

10 

20 

40 

60 

80 


Gojii utabional performance of edge-disQcxnt 
lienris cic~2 


o.lO 

0.15 

0,20 



0.01 


— 

0.02 

" 


0.03 

0.03 

0.04 

0.04 

0.04 

0.06 

0.13 

0.07 

0.09 

0.13 

0.17 

0.23 

0.44 

0.32 

0.81 

1.40 

0.4 4 

1.11 

2.41 


0.52 

1.18 

1.24 

1.86 

3.58 

3.35 

3.82 

6 . 50 

12,32 

3.20 

2.67 

3.07 

6.39 

9.83 

12.65 

11 . 14 

19.76 

34.02 


2.99 

3.71 

4.91 

lu 70 

17.95 

23.45 

2 u . 22 

28.55 

62.38 


0.30 

0.50 

0,80 

0.02 

0.03 

0.05 

0.03 

0.06 

0.12 

0.04 

0.09 

0.22 

C.IO 

0.18 

0.33 

0.21 

0.57 

1.29 

0.40 

1.00 

1.93 

c . So 

1.41 

2.55 

2.68 

5.06 

12.54 

5,64 

11.97 

27.44 

2.55 

4.64 


10.96 

20.03 


20.45 

44.14 

*■ 

4.44 



20.40 

... 


60 . 51 





- 


100 
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' i'aJale 5 . 3 . 3 : 

jjenBxt;/’ 

¥o. of nodes 

10 

2 u 

40 

60 

80 


C'^ rcacional perf orciajacL ol node-axs joint 

llo j-iofcic-1. 


t'.'iO 0.15 



0.03 

0 . 0 ^ 

0.04 

O .04 

0.07 


0.12 

0.14 

o.is 

0.71 

O '. 20 

0.72 


0 . 38 

2.12 

1.32 

4.23 

2 . 78 

7.71 


1.37 

4.82 

3 . oC 

18.70 

:i .62 

25.61 


3.10 

5.38 

24 . S 8 

54.61 

32 . 22 

56.64 


0.20 0.30 


0.01 

O.Gl 

0,02 

0.02 

0 .03 

0.03 

0.05 

O.OC 

0.06 

W .22 

0.12 

0,29 

0.55 

1.25 

1.51 

0 . o 7 

1.63 

3.01 

1.48 

3.60 

9.28 

23.80 

14.09 

31.60 

4.07 

14.37 

36.13 

65.14 

49.11 

94.17 

4.62 

93.23 

94,31 

- 


C.50 0.8U 


0.02 

0,06 

0.05 

0.08 

0.06 

0.12 


0.15 

0.67 

0.50 

1.54 

0.85 

1.94 

2.95 

7.33 

9.46 

21.80 

13.09 

28,05 

10.48 


43.83 


53.31 


- 



100 
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!i'a'':)le 5. 3. 4= Ooupuc tional perforijiance of node -ciis joint 
louris uic-2c 


ue. 3ity 

, of i'odes 

. .10 

0.15 

0.20 

0.30 

0.50 

0,80 

10 



0.01 

0.01 

0.02 

0.03 



- 

0.02 

t .02 

0.04 

0.06 




0.03 

0.03 

0.C5 

0.08 

20 


0.04 

0.04 

C . 06 

0.13 

0.16 


1. .03 

0.05 

0.09 

C.ll 

0.20 

0.31 


O.o7 

0.07 

0,10 

C.15 

0. ?6 

0.36 

40 

l .16 

0. 20 

0.25 

0.40 

C.65 

1 .07 


0.24 

0. 41 

0.53 

0.71 

1.03 

1.66 


0.27 

0.4S 

0.82 

0.96 

1.30 

2.14 

60 

'■.47 

0.85 

0.89 

1.56 

2.09 

MU 


1.05 

1.40 

1.55 

2.42 

5.1:> 

- 


1.33 

1.64 

2.^7 

3.09 

4.25 

*■ 

.'0 

e .90 

1.37 

2.04 

3.07 

__ 



?,1G 

3.13 

5.36 

3 . 00 

- 



.>.07 

3.86 

5.05 

11.18 



100 

-.33 

2,80 

3.94 





..99 

4.43 

8.78 


- 

- 


3 .36 

8.00 

7.99 
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NUMBER OB N0DES = i,0 ‘ 

NUMBER OF C0MM0DJT;:'5==5 

EDGE-DlSjOiNT HFURISTSC-' 
NODE-DISJOINT HEURISTIC- i 
EOGE-DIS JOINT HFURiSTiC- 
NODE-DISJOINT HEURISTIC- 


’ 0 


!'!G.b-3'1 COMPARISON OF HEURISTICS 


i 


I 



1x1 


OKj TfMP- ^ ^ .-r A T r^K x ^jrrQf^nc 


I 


NUMBER OP NODES ' 40 
'f hJ OF COMM ODl TIE S- 



VAR^Af 3^' OF CHMR TIME WITH NO- OF 
COMMtOiT.^-E 


r 



1 - .1 
os 0-8 1-0 


LFNSITV 

FiG 5‘3-3 VARIATION OF COMP. TIME WITH NO. OF NODES 



CHAPTER YL 


A"V"i>TlJ"JS POE POETIIBR EESEAt^ClI 

In oliifj -’■o’L an at’Lempc has been i''ade to propose 
soluoron teciPiineo lor a class of disjoinx flow problems, 
rlie proposed jao i, juS at'e not perfect, liorever they can be 
credited for no., ir Omg research in this field. Methods 
may yield bebLc.' ■ ocults by inbroduciny, soiie intelligent 
raouiJ’icationo. I'or example, a different upper bounding 
cr'iteria liiay bo sv, gesced for braxicn ano bound algorithms. 
Siaiilarly xiucn scoter performance can be expectea frora the 
comjjucer progra'-''- if they are suitably moaified. 

Jjevclopjient of exact methods for multicommodity 
ms joint nows shouH o be the concem of uhc persons 
interested in bins field. 

I'he preciei u \rork was concerned Tribii jiure Disjoint 
Llox^, i.e., wliure oJ-l ai’cs or nodes allow only one commoaity 
to flow over it. It inay be wortlruhile to coiioidei the inxed 
aisjoinfc fl oxrj, i.o , where some arcs or nodes allow 
ci"ialtaneous Ilou oi commodities. 
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APPEEDIX A 


A f azid mol/liods y ‘wiiicii s,r0 used in ^116 

V( ■. . 4 , -’i , 

algorithms 'd‘f 'cah .■jilo'\?lous , chapters ,‘'*"are hescriled here. 
■A-1 'SIH(JI,E _q0Ii^:0*;Ii'Y MAZMAD i!L0¥ ilG-Ohll’m 


'ifinds maximal iloi-r fi'om a soui’ce node 

'-Saah node iia capaoitated network. Algorithm 

' ■ ' ■ '■ ■ ' ' ^ ^ V‘ 

proposed hy Pord and Pulkerson [8] is used in this work and 

described here briefly. Por other algorithms for the same 
problem see [4,21]. 

ii'lie algorithm begins with an initial 





I': v ;, i 


vector Z. Some loioii'm feasible flow vector or X=0 cj 


used as che inr uial flow vector, let p be 
feasible flow vector 
be divided into two routines; labelling routine 


improvement routxne. labelling routine is' a sysremaTic 


method for che elding whether a flow augmenting- path with 


respect vo X exists or not. Plow impi'ovement 


che routine to improve the flow on the augmenting path 
discovered by t]ie labelling routinb.- 

ksbdl lldg Routine; : /.-siabel the source' an (s'*',00) and 

^ ‘ , ' ' , ‘ U ‘ I' f. 

continue labolliag thf-fo^^^ion^^laboll^ rules; 






(1) If 1 is a labelled node and 3 is an unlabelled node 

at this staj^e and (i,i) e A and ^ label j 

with (i , e .) where e . - min [e. , c. . - x. .]. 

y y J- , -X. J 

(2) If i is a labelled node and j is an imlablled node’ 

at this.Bta(3,e ahd (3,i') e A and - y 0, label j 

with (i”, e,) where e. = rain fe. , z.. I- 
J ' 3 '-1 ji-’ 

•i ; -t--' ,.f , ■ , . ’ 

I'low Improveme nt Soutine : Assume that sink is labelled 
{Jc e^) 'Ee place by x-^ + e^ and turn to 

If Ejj. is labelled as ( 3”, e^) , replace z^^. by z^^^.-e^ and 
turn to N.. Contitiue until source is reached. Erase all 

J 

the labels and start labelling routine. 


labelled nodes are scanned to label the unlabelled 
nodes in the oiuer they arc labelled. Iflien sink node is 
labelled, the labelling routine terminates and flow improve - 
riient routine .be giiifs. If all the labelled nodes are 
scanned and sink is still unlablled, the algorithm stops. 

Ihe present flo-fr is the optimal flow. - • , 

A, 2 EAXCMAI ACfCIiOHlINC} PilH ALGORITHM , ^ 

Maximal augmenting path in a network is' ^ tbo augmenting 
path on which ma^'iimum amount of flow can be' .augmented. The 
flow value that can be augmented on a path is laioim as the 
augmenting capacity of that path. The method to find such 


a path io simiXar to 




i)ijkstra*s method [ 5 ] to solve shortest 


path 








m 


" ' K> 


Classify the source node as node xjith permanent label 
and assi^ it a lal^cl +O0- All other nodes have temporary 
label -OQ. lei acurce node be the peinianently labelled 
nods iidd©i*-,looi<»id.#ratt!ion. 


Step^l Chance the tentative labels p-i’s of the temporai’ily 
labelled nodes 3 employing the fol3-Oin.ng scheme, xfhere i 
is the permanently labelled node mder consideration 
with label u. : 


mn 


min 


St ep_.,2 % determine the largest of all tentative la 

and classify tliat node as permanently labelled' node 

■ v..‘, '-.r 

this node now be the node under consideration. ''Qo' 


c; jLo 




When, after at ao^t 
tal iterative stejs^tsin: 
procettur© t©minat|^,^. ,..,1 










83 







A. 3 BRANCH JHIi) BOUNi; APPROACH 


Brancli and boimd approacH is a very useful tool for 
solving combinatorial optimization problems. It provides 
a systematic method to search for the optimal feasible 

’tf j_ >' 

solution by doing only partial enumeration. In the coursb 
of applying the branch and bound approach the overall set 

partitioned into many simpler 
one promising subset is chosen and':- 
;an>ffort to- find the best feasible solution from;;;;! 

^hen^ that subset is said to be fathoid, 
^thelbeli'^lasible solution is not _ found then that subset 
is again paruirioned into' two or more simpler subsets 
(this opera'oioii knowa. as branching) and the same process 
is repeated again. At each partitioning of a set’ into * 

i'-' j , ' ' f ff ‘ 

ibs subsets a lower bound (or upper bound) , which gives' the 
niniffium value (or maximum value) of the objective function 
/rhich any -feasible solution in that subset may yield. 

d,, , ''•s'/ ' ' 

this lower bound is greater than (or lower than) some 
best knovm feasible solution (known as incumbent) the: 
that subset is not further branched (tliis operation to 
as pruning), Ifnon all subsets are either fathomed or 
pruned, incumbent provides the value of an optima.! ^ o': 













<lil|l|i«, N00gji;o , OEMSf If , ARCNO , SOURSE , SI NK , START , 


€♦♦♦** 'THIS WRt CeWERATES THE SODE^NOOi INC I01WCE'; MATRIX 
DO 10 IisiVKODBNO - ^ ' ,’‘,. 4 , 


♦♦♦ #^#111 )jE f.)({ f f f !^ ♦ f ;j! f f ♦** ijt 4 : f ♦ t 3 ^ Jliffl: ) 

%Ni|jif«$"A CONNECTED NETWORK 


tM t N0N8ER HE NODES IN f«E NETWORK 
i$|f|ER; s nOMBSR or COMMODITIES EI^OWING 
.«l|»*ls|S|III»»EltfOr ARCS IN THE NETWORK 

DENifX { A t.OfiHri? WHICH CONfROM THE DENSITY OP THE NETWORK 

SODRSE(i) ; Si.-RR^;fc, NODE OP COMMODITY I 

SINKtr); SI-''!* f-ODt' OP COMMODITY I ; ' ' 

STURfU) S STARTING NODE OF ARC 0 

EWDCaa i ENDING NODE OF ARC J 

WU-)h CAFACin Of ARC J ' .r" 

^|ljNeCf I THE NODE^NQDE INCIDENCE MATBIX GENERATED BY THE PROGRAM 

iTNTEOiR CONECTdOO, 1003, DlSfC too), START<2000),iNDC2000)r 
1 CAP C2000 3 rSOURSeC 10 ) , SIHKC 1 05 , RANDl ,RAHa2 , TEMPI , ARCNO 

;.v:; 

. ■* '' ,•' ‘ ' > 

' ■■— ■.■■•-r-p . ■ - 


DO ^'S 4*sl,NODEHO 
S CONiXff (i*a)«o 

1!^ MST(I5«0 • 

t>0 15 t«3,NO0ENO 
ir(RAN{X).0f.DENSTY)GO TO 15 
!.,=l,+ l 
f’or.,i iKUE 
l>u ^0 0 = 2,1, 

COliLC'i f i ,J)=1 

Cbr.rc I U! , i )=-l '. 

L 3 ^'’rf'’b=ij 3 STNDi- 1 
lasi (.j i = i 
r b I f T i< i : r; 

}Hi vs ]=2,NUlitliO ; 

FUl'K = u . ■ ' 

UU V5 J = 1 , MUOffUJ i- 

JT (v),k,(3.3 .(!H,Cnj^PCTCt,J3,NE,0)GO TO 25' 

Tf f Kouc.fM.nGo Tu io 

1 : r ( i. J S T H ( 1 , h y , N 0 u E « 0 . 0 H , D I S T ( J ) , E Q , I ) GO ’ : TO ■ I # ; ■ ' 

i,.»r»‘iNC}sLO sTM'j+1 . ■ ' 

I .Gl'.UKwSTlilGU xn 25 





mim. 


iff *OiJi3fe.RikND2, £0,0)60 to 40 

” tii$f <Rji«o t )«iifsT f ) 

. .I<i$i(RA'.'i02J«f£«P1 ' ■ 

■ .. , ,, , . ' ■ ' ' 

^.•■■•f'lilS PART' FIMDS P0SITI0H5 OF S0UR5ES AND SINKS 

I it NUMBER 

ifir ^-ppRSEClJtsRAMajtNODENO 
;'St»KCJ3»RANCX3'»NO0Er»O 

IfCSOyRSECI).EO.O'«DR.SINKa).EQtO.OPtSOiJRi5E(I>,EO,5INKCl))GO TO 52 
' 'PNTINUE ^ 

m**'** ARCS ARE lOENTIFlEO FROM THE INCIDENCE MATRIX 

DO SO I»1,P»EN0 " . ■ •.-■ ■ ' 

DO §0 J»ll,NODENO 

■ IFtCONECT-C 1 ,iSi .06,0 >0Q. TO, 10 , ,, 

‘ ^ 

START CK 1 «I.IST C I ) ' %>■. ■ 

BROCK 5»tisTcj3 ' . 

AO CONTINUE ‘ ^ ' 

ARCNOaK 

GENERATE CAPACITIES FOR ARCS . . 

tiO 70 Jsl,ARCt«0 , : ;• 


C APtJ)sttRANCXJ*lOO} 
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